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Abstract. The B a meson is an interesting particle to study because a sizable mixing induced CP 
violation in the B s — B s system would be an indication for physics beyond the Standard Model. In 
this paper we present a measurement of the lifetime difference AF between the B a mass eigenstates 
and the CP violating phase in the decay B s — > J/tp<j>. In 1.7 fb of data collected with the CDF II 
detector at the Tevatron pp collider we measure Ar = 0.076+°'^ (stat.) ±0.006 (syst.) ps" 1 , well 
consistent with the Standard Model prediction, and a mean B 3 lifetime of ct s = 456 ± 13 (stat.) 
±7 (syst.) /im. We find no evidence for CP violation [U12]. 

PACS. 13. 20. He Decays of bottom mesons - 14.40.Nd Bottom mesons 



1 Introduction 

In the B s -B s meson system the flavor eigenstates are 
not the same as the mass eigenstates. The mass dif- 
ference between the heavy and light mass eigenstate, 
B s h and B s l, determines the frequency of the oscil- 
lation of the B s mesons. Two other quantities which 
affect the time evolution of B s mesons are the decay 
rates I# and -TV of the two mass eigenstates. The dif- 
ference Ar = /Y — I# was measured first by CDF [3] 
and recently with higher precision by D0[4]. 

If the difference Ar is larger than a few percent of 
the mean decay rate r = {Tl + Ih)/2 a time depen- 
dent analysis of B s decays without flavor tagging be- 
comes sensitive to a further quantity, the CP violating 
phase (j)g. This phase describes the mixing induced CP 
violation and is related to the angle f3 s in the nearly 
degenerated unitarity triangle obtained from the mul- 
tiplication of the second and third column of the CKM 
matrix. The Standard Model expectation value for <p s 
is very small [5 . Therefore a measurement of the phase 
which deviates significantly from zero would indicate 
new physics. 

To determine Ar the lifetime distribution of B s 
decays is measured. Because it is very challenging to 
distinguish the two components of the lifetime distri- 
bution additional information is needed to separate 
the light and heavy mass eigenstates. Therefore we 
exploit the fact that the two mass eigenstates are re- 
lated to the CP eigenstates. In case of no CP violation 
{<j) s = 0) B s h is CP odd and B s l is CP even. 

A decay mode that allows to measure both lifetimes 
is B s — > J/4>4> with J ftp — > and <p> — > K + K~ 

which is a composition of CP even and odd states. 
Because the B s is a pseudo scalar and J/ip and <p are 
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Fig. 1. Definition of transversity angles 0, <p and tp. 

vector mesons, the orbital angular momentum between 
the two decay products can have the values 0, 1 or 2. 
S- and D-wave decays are CP even, P-wave decays are 
CP odd. Consequently, the two CP eigenstates can 
be separated by their different angular distributions 
of the decay products. 

The angles u> = (cos 6, (p, cos ip) used in this anal- 
ysis are defined in the transversity basis illustrated in 
Figure [TJ and <p are the polar and azimuthal angle of 
the [i + in the rest frame of the J ftp where the ir-axis 
is defined by the direction of the B s and the xy-plane 
by the cp — ► K + K~ decay plane, ip is the helicity an- 
gle of the K + in the cp rest frame with respect to the 
negative B s flight direction. 

2 Data Sample and Selection 

The analyzed data sample with an integrated luminos- 
ity of 1.7 fb -1 was collected by the CDF II detector at 
the Tevatron which collides pp at a centre of mass en- 
ergy of 1.96 TeV. The detector components essential 
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Fig. 2. Network output for B s signal and background 
events. 



for this analysis are the silicon vertex detector for a 
precise lifetime measurement, the central drift cham- 
ber for good momentum and mass resolution and the 
central muon chambers for the identification and se- 
lection of muons. In addition to the energy loss in the 
tracker the time of flight detector is used for particle 
identification. The B s — > J/ip<fi events are triggered 
by a pair of two oppositely charged tracks matched to 
signals in the muon chambers and having an invariant 
mass close to the J mass. 

The J ftp candidates are combined with 6 candi- 
dates built from pairs of oppositely charged tracks as- 
sumed to be kaons. The B s candidate is obtained in 
a common vertex fit. After applying basic kinematic 
cuts a neural network is used to improve the selection. 
The network is trained on simulated B s signal events 
and background events from B s mass sideband data. 
Kinematic, particle identification and vertex fit qual- 
ity variables are used as input to the network. Figure 
[2] illustrates the good separation of signal and combi- 
natorial background events. 

With a cut on the network output that optimizes 
the signal significance 2500 B s — ► J/ip(f> decays are 
selected. 



is given by 

d i P(u,t) 
dwdt 



c< |A)l7iMT + + |^||72MT + 

+ |^| 2 / 3 (^)T_ + |4)||A|||/ 5 Mcos(<J||)T + 

+ |A||||A±|/ 4 McOB(4l-*||) 



sin <p s (e' 



-r H t 



-r L t 



)/2 



\A Q \\A x \f 6 (u)cos(S ± ) 



sin 4> s {e 



-r H t 



,-r L t 



)/2 



(1) 



with 



,-r L t 



2cos 2 ^(l - sin 2 6 cos 2 



sin 2 9 sin 2 . 



((l±cos 

9 

32^ 
9 

sin 2 tp(l 
— — sin 2 ib sin 2 

327T 

9. 2 . 
— — — sin ip sin 29 sin 6 
32n 

9 1 2 

= sin 2ib sin 9 sin 26 

32ir V2 

9 1 

= sin 2ib sin 29 cos 6. 

327^ 



+ (lTcos0 s )e-^*)/2 



Note that this distribution is invariant under the trans- 
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Because of this four fold ambiguity this measurement 
is insensitive to the sign of both, 6 S and Ar . 

The finite lifetime resolution and differences of it 
between signal and background are included in the fit 
model. The angle dependent acceptance is taken into 
account by an acceptance function obtained from sim- 
ulated events. The good description of data by the 
simulation is exemplarily shown for the selection net- 
work output in Figure [3J 

The fit projections for mass, lifetime and the angle 
cos tp are shown in Figure [U 



4 Result Assuming no CP Violation 



3 Mass, Lifetime and Angle Fit 

The mean B s lifetime cr s = c/T, the lifetime differ- 
ence Ar and the decay amplitudes A)> Aii , and A± 
of the three angular components with relative phases 
Sj_ and 8i\ are extracted in a 5-dimensional unbinned 
maximum likelihood fit in mass, lifetime and angular 
space. Empirical models are used for the background 
distribution. The B s mass signal is described by a sum 
of two Gaussians. The lifetime and angle distribution 



We first consider the case of no CP violation (6 S = 0). 
This simplifies the fit model because the last two terms 
in equation ([1]) and the fit parameter Sj_ vanish. 

The result can be affected by several systematic 
uncertainties which are evaluated using pseudo exper- 
iments. The investigated effects are the influence of 
the angular background, the signal mass and the life- 
time resolution model, the contamination of misrecon- 
structed B° — > J/ipK* decays, the acceptance function 
and the silicon detector alignment. The largest system- 
atic uncertainty for Ar is the B° cross feed and for 
cr s the lifetime resolution model and the alignment. 
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Fig. 3. Distribution of the selection network output for 
sideband subtracted data and simulation. 



Setting <p s — in the fit we obtain 



CT S 

Ar 

|2 



456 ± 13 ±7 urn 
0.0761^63 ± °- 006 P s ~ 



|A)| = 0.530 ±0.021 ±0.007 
|A||| 2 = 0.230 ± 0.027 ±0.009 

The first is the statistical and the second one the sys- 
tematic uncertainty. Since the likelihood scan for the 
strong phase 5\ \ , shown in Figurc[5l has a non-parabolic 
shape due to a symmetry at (5|| = tt we do not quote 
a point estimate for this quantity. 



5 Fit with Floating CP Violating Phase <p s 

Since a maximum likelihood fit is only guaranteed to 
be unbiased in case of unlimited statistics we studied 
the fit in pseudo experiments. In case of free (f> s param- 
eter we observe that for low input values of Ar or 4> s 
there is a bias towards higher values. This is illustrated 
in Figure [U 

The bias can be understood by looking at equation 
(fTJ) . If 4> s approaches zero the last two terms vanish and 
5± becomes undetermined. This means that the fit can 
not improve the description of the data any more by 
varying <f> s . It effectively lost a degree of freedom. The 
situation is similar when Ar gets zero. Then <f> s and 
5± are undetermined. 

Because of the biased fit result we do not quote a 
point estimate for Ar and <p s , but construct a con- 
fidence region following the procedure suggested by 
Feldman and Cousins Ej. 

For each pair on a grid of assumed true values of 
Ar and <j) s we calculate a p- value, which quantifies the 
probability to get the fit result observed in data. To 
determine the p-value we use the likelihood ratio 
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Fig. 4. Fit projections for mass, lifetime and angle cosip. 
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Fig. 5. Likelihood scan for the strong phase 5\ 
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Fig. 6. Distribution of fitted AT (left) and cj> s (right) in 
pseudo experiments. The blue arrow indicates the input 
value. 



where 9 are the nuisance parameters and the hat in- 
dicates the parameter values which maximize the like- 
lihood C The R distribution for assumed true values 
of AT and </> s is obtained from pseudo experiments. 
The pseudo experiment input values of the nuisance 
parameters are taken from a fit to data, a procedure 
known as plug-in method. 

The p- value is the fraction of pseudo experiments 
with an R value higher than the one in data. The 90% 
(95%) confidence region is then defined by the Ar-cf> s 
points with a p- value above 10% (5%). 

The result is presented in Figure El Note that only 
the first quadrant is shown. The other three quadrants 
can be obtained via the transformations given in equa- 
tion ©. The p- value for Ar = 0.1 ps" 1 and <j> s = 0, 
which are approximately the Standard Model predic- 
tions, is 22%. 



6 Conclusions 

Using 2500 B s — > J/ip(f> decays selected by a neural 
network in a data sample of 1.7 fb -1 CDF has per- 
formed a mass-lifetime-angle fit to measure the lifetime 
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Fig. 7. Confidence region in the Ar-<f> 3 plane. Only one 
out of four possible solutions is shown. The others can be 
obtained by the transformations <j> s — > — <j> s and Ar — » 

—Ar, A s — > 4> s + 7r 



difference Ar between the B s mass eigenstates. The 
value obtained under the assumption of no CP viola- 
tion is consistent with the Standard Model expectation [5j 
and previous measurements 3,4]. The extracted mean 
B s lifetime is currently the most precise measurement 
and agrees well with the world average B° lifetime as 
predicted by theory. 

If we allow for CP violation in the fit model we 
observe a bias away from low Ar and <j) s values. It is 
understood by the structure of the likelihood function 
and limited statistics. Instead of a point estimate a 
confidence region is determined in a frequentist way. 
The result is compatible with the Standard Model 
and can not rule out any minimal flavor violating new 
physics scenario which changes the phase <j> s , but does 
not significantly affect b — > ccs tree level dominated 
processes. 
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